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Abstract 



high-frequency oscillations of Newton's constant, G. The effect occurs because 



inflation can shift the expectation value of a non-minimally coupled, Brans- 
Dicke-like field away from the minimum of its effective potential. At some 
time after inflation ends, the field begins to oscillate, resulting in periodic 
variations in G. We find conditions for which the oscillation energy would 
be sufficient to close the universe, consistent with all known constraints from 
cosmology and local tests of general relativity. 
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I. INTRODUCTION AND SUMMARY 



Time-variation of Newton's constant, G, has been the subject of many investigations. 
Most studies, motivated by Dirac's large-number hypothesis, focus on monotonic variations. 
Brans-Dicke theory, for example, was devised as a field theoretic model that incorporates 
a monotonic variation in G. Its key element is a massless, scalar field, 4>bd, which couples 
non-minimally to the scalar curvature, R, through the interaction, <PbdR- The effective 
Newtonian constant is proportional to <\>b\>- The massless 4>bd increases uniformly in a 
matter dominated epoch, resulting in a monotonic decrease in G. The model is problematic 
from a quantum field theory point of view because there is no reason why the mass of 4>bd 
should remain zero after quantum corrections. The problem from an empirical point of view 
is that there are tight limits from cosmology and local tests of general relativity (for a review 
see i) 

In this paper, we consider the possibility of an oscillating Newtonian constant in which G 
varies with a frequency much greater than the expansion rate of the Universe, v 3> 10~ 17 Hz 
0. A model with oscillations in G is quite simple to construct; e.g., a massive Brans- 
Dicke model in which a scalar field $ (related to <Pbd) is displaced from the minimum 
of its effective potential. The displacement in <3> could result naturally from an epoch of 
inflationary expansion. During inflation, the scalar curvature is non-zero, and so produces 
a substantial contribution to the effective potential for $ through the $i? interaction. The 
added contribution shifts the minimum of the effective potential for $. After reheating, the 
scalar curvature is zero, and <3> is free to oscillate about the true minimum of its effective 
potential. From the point of view of quantum field theory, the massive Brans-Dicke field 
does not suffer the naturalness problem of the massless case. One might anticipate that the 
mass of $ would be similar to that of typical elementary particles, m > 1 GeV or so, a range 
corresponding to very high-frequency {y~ l ~ mT x << 1 sec) oscillations. 

Accetta and one of us have previously discussed some effects that high-frequency 
oscillations of G could have on cosmological measurements. The purpose of the present 
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study is to examine a very plausible scenario in which inflation sets off the oscillations in 
G. For simplicity, we consider cosmological models for which the energy density contains 
three components: ordinary, baryonic matter, p#; radiation p^, and oscillatory energy, p$. 
It is straightforward to extend our results to models which include a fourth, conventional 
dark matter component. We assume that some inflaton field completely independent of $ is 
responsible for providing the vacuum energy density needed to drive inflation. We further 
assume that $ is so weakly coupled that ^-particles are not produced by the decay of the 
inflaton during reheating. During inflation, the $ field is shifted from the minimum of its 
effective potential; after inflation, the only effect on oscillations is the Hubble red shift of the 
$-field kinetic energy. We find that the scenario passes all known tests for a wide range of 
model parameters. In particular, it is possible to find viable models in which the oscillatory 
energy in <3> provides all of the missing energy density needed to reach the critical density 
without introducing conventional dark matter. 

The organization of the paper is as follows: in Section II, notation and the basic equations 
are introduced to describe massive Brans-Dicke cosmology with oscillating G. In Sections 
III and IV, the oscillations in G induced by an epoch of inflationary expansion are computed 
for scalar fields coupled to the scalar curvature through linear and quadratic interactions, 
respectively. The red shifting of the oscillation amplitude after inflation is analyzed to 
determine the oscillation energy today. It is shown that the oscillation energy is compatible 
with all cosmological constraints for a wide range of parameters and that it is possible for the 
oscillation energy to reach closure density. In Section V, the constraints from local tests of 
general relativity are examined, including measurements of light deflection and time delay, 
data from laser ranging to Earth-orbiting satellites and the Moon, and tests of Newtonian 
gravity using tower and laboratory experiments. Section VI presents concluding remarks. 

II. BASIC EQUATIONS 
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A. Action and field equations 



The gravitational theory underlying the oscillating-G cosmologies is a general scalar- 
tensor metric gravitational theory, in which the scalar field has both a mass and a non-linear 
A0 4 self-interaction. The action has the general form 

I = J ^rf 4 x[M 2 /(0)^-^V / ,0V>-^(0)] + / m (^ m ,^), (2.1) 

where M 2 = (167rG) _1 , g^ is the physical metric, I m is the matter action, with matter fields 
ip m coupling only to g^v (hence a metric theory) , and R is the Ricci scalar constructed from 
g^y. We use the standard notation and conventions of £j],||]. The potential V is given by 



V{<f>) = \m 2 <? + ^A0 4 . (2.2) 

With V = 0, this would correspond to a generalized Brans-Dicke (BD) theory |l] with 
cj) BD = /(</>) and the BD coupling function 

u> = f/2(Mf') 2 , (2.3) 

where "prime" denotes d/dip. The field equations are 

1 ~ 1 1 1 

G> = 7^pJ T ^ + j(f;^ ~ 9^ n /) + 2M2j(M," - 2^,a0' A - 9^V((j))), (2.4a) 

« , g ~ 2/V - jf(l + 6M 2 r)0, A 0- A + Iff 

" / + 3M 2 (f) 2 ' 1 J 

where is the physical stress-energy tensor derived from the matter action I m (ip m , g^), 
and □ is the d'Alembertian using the physical metric, and so on. 

For cosmology, it turns out to be useful to recast the theory into a representation with 
an Einstein gravitational action by making a conformal transformation to the non-physical 
metric 

9tiv = f{4>)g^- (2.5) 
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This representation is sometimes called the Einstein conformal frame (see for discussion). 
The action then takes the form 

I = J V^d 4 x[M 2 R - iv„ W$ - V($)] + I m (ipm, /"V). ( 2 - 6 ) 

where V(<&) = f~ 2 V((j)) and $ is defined by the differential equation 

d^/dcj) = r\f + 3M 2 /' 2 ) 1/2 - (2.7) 

Damour and colleagues [^-[7j have discussed general classes of scalar-tensor theories in this 
framework; their scalar field <p is given by <p = $/2M, A(<p) = f~ 1/2 , and V = 0. We 
will consider models in which the non- minimal coupling is linear, /(</>) = 1 + a>4>/M, and 
quadratic, /(0) = 1 + £<p 2 /M 2 . 

The field equations in this representation are 

G> = ^27^ + A - ^m^,a$* A - <^(*))> (2.8a) 

□$ = ^(n$)-^r 2 T). (2.8b) 

where now □ and G^ v are defined using the non-physical metric; however, T^ v still refers to 
the physical stress-energy tensor defined using locally measured quantities and the physical 
metric. 

B. Cosmology 

Assuming a standard homogeneous and isotropic universe, we choose coordinates (t, x) 
so that the "line element" constructed from has the Robertson- Walker form 

ds 2 = -dt 2 + a(t) 2 d(r 2 , (2.9) 

where a(t) is the scale factor, and da 2 is the spatial line element representing the hypersur- 
faces of homogeneity. Since we will assume a standard inflationary scenario for the early 
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universe, we will adopt a spatially flat metric, i.e. we will choose k = 0. The relation 
between t and a and the corresponding physical variables t and a in the Robertson- Walker 
version of the physical metric g^ u is dt = f~ l l 2 dt, a = f~ l ' 2 a. The physical density and 
pressure p and p have their usual meanings as local properties of the matter measured in 
local Lorentz frames of the physical metric g^ u . Defining H = a/a, where an overdot de- 
notes derivative with respect to t, we obtain the cosmological field equations and equations 
of motion 

H 2 = (6M 2 )- 1 (^pf- 2 + U> 2 + V(<$>)) , (2.10a) 
£ + 3H$ = -dV{<S>)/d<& - i(p - 3p)d(f- 2 )/d<$>, (2.10b) 

p=-3(p + p)(H-±f/f). (2.10c) 

Note that the use of the non-physical metric g is the origin of the pf~ 2 and the —\f/f terms 
in Eqs. (|2.10a|) and (|2.10c|) , respectively. The physically measured Hubble parameter H = 
ar x ddjdt is related to the Hubble parameter used in the Einstein conformal representation 
by H = f 1 / 2 (H -///). 

If p — 3p is ^-independent, as we will assume in this paper, we can use Eq. ( |2.10b| ) to 
define an effective potential 

v cS = v^) + l -{ P -3p)r 2 
= r 2 (v(<p) + \(p-sp)), (2.ii) 

such that $ + 3if$ = -dV cS /d§. 

We assume that inflation is triggered by some inflaton field in the matter sector (not $!), 
whose vacuum energy density pv dominates the total energy density and leads to an epoch 
during which py = —p = constant. We relate the vacuum energy density to a temperature 
by the standard relation (in units where h — c — ks — 1) 

p v = T V = (10 14 GeV) 4 ^ (2.12) 
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where t v = T y /10 14 GeV. Note that, with these definitions, p v /6M 4 = 2 x 10^%. 

In characterizing cosmological evolutions, it will be useful to define a number of key 
variables. 

(a) Oscillation amplitude: 

!a<j)/M, [linear] 
(2.13) 
i<p 2 /M 2 , [quadratic] 

(b) Coupling parameter: 

!3a 2 , [linear] 
(2.14) 
12£ 2 2 /M 2 , [quadratic] 

(c) Anharmonicity parameter: 

C = A0 2 /2m 2 , (2.15) 

(d) Ratio of scalar-to-radiation energy density 

x = p^/r 2 p R , (2.16) 

where pn is the energy density of radiation, and 

P$ = U 2 + V($) . (2.17) 

The parameter ( determines the relative importance of the A0 4 anharmonic potential com- 
pared to the mass term m 2 (f> 2 . Note that (5 is related to the Brans-Dicke parameter uj by 
uj = |(1 + rj)/P; as (3 — > for finite rj, oo — > oo, and the theory becomes effectively general 
relativistic. 

For future use we also define the dimensionless coefficients 

(e) Ratio of Planck mass to scalar mass: 

y = M/m = 1.7 x 10 18 (lGeV/m) , (2.18) 

(f) Q 

Q = \2\-\p v /QM A ) = 2.4 x 10- 17 A _1 4 . (2.19) 



We assume that the material content of the universe is strictly that of ordinary matter (no 
exotic dark matter), and we will treat the equation of state for matter as having three phases, 
inflationary, radiation dominated and matter dominated, with the ratio p/p = v = constant 
during each phase, with the respective values -1, 1/3, and 0. 

In all our models we will assume that the present amplitude of G-oscillations is very small, 
i.e. t]q « 1, where the sub- or superscript denotes present values. We define density 
parameters at the present epoch in the Einstein conformal frame: for matter (baryons 
plus radiation): f2 m ° = p° /QM 2 Hq; for radiation alone: £l R ° = p° r /QM 2 Hq; for the scalar 
field: = p$° /6M 2 Hq. Because of the present smallness of 77, and because of its rapid 
oscillations compared to the expansion timescale, the physically measured H ~ H , so fi^, 
fl° R , fl^, and Xo are interchangable with the physically measured quantities. We note that 

Q° R = 4.3 x l(T 5 /i~ 2 (T /2.73) 4 , (2.20) 

where h is the present Hubble parameter in units of 100 km s" 1 Mpc -1 (0.4 < h < 1), and 
T is the temperature of the cosmic microwave background in degrees Kelvin. 

Since $ interacts only gravitationally , it is thermally decoupled from ordinary matter; 
we ignore the possibility of particle creation induced by the non-minimal coupling of to 
curvature. Consequently, we can relate the energy densities of baryons ps and radiation 
p R to the scale factor a and temperature T using the matter equation of motion in the 
physical representation, dp/dt = —3(p + p)H, combined with standard thermal equilibrium 
arguments. Thus, when baryons and radiation are decoupled, we have p R oc a~ 3 , p R oc a -4 , 
T R oc cT 1 . Then, denoting the end of inflation, marked by the conversion of false vacuum 
energy density to radiation, by I, it is straightforward to show that the non-physical scale 
factors at these epochs are related by 

oo/o/ = (T / /T )(/ // / ) 1 / 2 « 4 x 10 26 (/ ///) 1/2 (Ty/10 14 GeV) . (2.21) 
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C. Evolution of oscillating scalar field 



In the limit of small the basic oscillation frequency of $ is ~ m. Roughly speaking, 
in an expanding cosmology, $ will oscillate when the period of oscillation ~ m _1 is small 
compared to the expansion timescale ~ H" 1 , i.e. when m >> H. Because the equation for 
$ is non-linear, the effective frequency of oscillation could differ substantially from m. As a 
consequence, we will find it necessary to give a refined criterion for the onset of oscillations, 
separately applicable for the linear and quadratic models which we are considering. Once 
oscillations commence, however, the subsequent evolution of the oscillation amplitude can 
be found using a generalization of the method used by Turner |§ . 

We first convert the field equation for $ into an approximate "energy conservation" 



equation, by multiplying Eq. (|2.10b|) by / n $ and extracting total time derivatives. The 
result is 

j t [f n {\& + V($))] = -3Hf n $ 2 + 6M 2 H 2 df n /dt, (2.22) 

where n = |(1 — 3z/) = 2,0,4, for inflation, radiation domination and matter domination, 
respectively, and where we assume that n is constant during each phase dominated by the 
corresponding form of energy. We define 

E = f n (h 2 + Vm, (2.23a) 

, . (2.23b) 

y/2(f-»E - Vm 

Then for periods << H , we average the right-hand-side of Eq. (|2.22|) over a period, 
with the result 

(3H(f n $) 2 ) ~ (3H/T) £ f n &dt = (3H/T) j> / n $rf$, (2.24a) 

(6M 2 H 2 df n /dt) ~ 6M 2 i? 2 T -1 [ T (df n /dt)dt ~ 0. (2.24b) 

Jo 
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Defining 7 = (1/ET) § f n $d$, we find 

E^dE/dt = -37# = -3jd/a . (2.25) 

For 7 ~ constant, this equation can be integrated to give E oc a~ 37 . Defining the 
integrals 

I ± = j>(\ - f n V(<5>)/E) ±1/2 f n/2 d$, (2.26) 

then 

7 = 2/ + /7_. (2.27) 

It is useful to write I± in terms of 0: 

i± = /(i - r-V(0)/E) ±i / 2 (/ + / 5) i / 2 r /2 - 1 #. (2.28) 

Note that when $ reaches its extremal values within each oscillation, $ = 0, so E = V(§ m ) = 
/( < / > m)~ 2 ^(0m) ) where the subscript m denotes the extremal value. In the quadratic coupling 
case, there is a symmetry — > — 0, so the maximum and minimum values of are identical 
apart from a sign; in the linear coupling case, this is no longer necessarily true. 

However, for oscillations generated by inflation, we will find that 77 << 1 already at 
the onset of oscillations, so that in evaluating 7, we can set / ~ 1. We can now treat 
<t>max — —4>min = 0m m ah cases. We define Cm = C(4>m), and x = 0/0 m . We also have, in 
the linear case, (3 = 3a 2 and in the quadratic case, (3 m = /3(0 m ). Then, independently of n, 

J± linea r = ^ + g^l/2 f 1 ^ _ ^±1/2, j + ? q ^ 

JO 1 + Cm 

^quadratic = 4 f X Q _ ^±1/2^ + _Cm 2)±1/2 (1 + ^2)1/2^ (2>2gb) 
JO 1 + Cm 

Table | shows the resulting values of 7 for the appropriate ranges of Cm and [3 m . For the 
linear case and the quadratic case with small coupling parameter f3 mi E has the expected 
variation: as matter (oc a~ 3 ) in the harmonic regime, and as radiation (oc a -4 ) in the 
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anharmonic regime. However, for the quadratic case with large coupling parameter, E 
for anharmonic oscillations varies like matter (oc a~ 3 ), while for harmonic oscillations, it 
decreases more slowly than matter (oc a~ 2 ). The evolution of the amplitude <j) m can be 
obtained from the fact that E « |m 2 m 2 (l + Cm)(l + Vm) 11 ^ 2 - Table | shows the relevant 
values of S, defined by m oc a~ 3S . 



III. LINEARLY COUPLED OSCILLATING G 

A. Oscillations generated by inflation 

Inflation can naturally lead to a finite displacement of $ from zero, leading to finite 
oscillations. We assume that a false vacuum epoch occurs during which — p = p = py ~ 
constant. Then $ + 3if$ = — dVi($>)/d$>, where the effective inflationary scalar potential 
Vi is given by (Eq. ( flUD ) 

A } (l + «0/M) 2 { ' 

During the inflationary epoch, the effective damping caused by the 3if $ term leaves the 
scalar field at the minimum of Vj, which occurs at dVj/d^ = = dVj/d(p. 

Figure [I] shows schematically how the effective potential during inflation, V/($), differs 
from the effective potential today. The effect of the nonminimal couplings is to change 
the shape of the potential, shift the minimum away from $ = 0, and raise the value of 
the energy density at the minimum. The change in shape depends on the nature of the 
nonminimal coupling; e.g., compare the upper and lower panels. For the linear case, we can 
use Eqs. fl2.18|) and (|2.19|) , to obtain the equation for the minimum, expressed in terms of 
the amplitude of the oscillations in /, 77 = a<ft/M, 

m {a 2 /\y 2 + r 1 2 I + l -r l ])=a'Q. (3.2) 

The first term in Eq. ( pj.2[ ) comes from m 2 in dVj/d(f), the second comes from A0 3 , while 
the third comes from a cross term between \<j) 3 and axf)/M in the denominator of Eq. ( [3.1|) . 
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Note that, in this case, the anharmonicity parameter £ (Eq. (|2.15|) ), which is essentially the 



ratio between the second and first terms on the left-hand-side of Eq. (|3.2j), can be written 



One class of solutions has the first term in parentheses in Eq. ( |3.2| ) dominant; this cor- 
responds to (j << 1, i.e. to scalar fields in the harmonic regime during inflation. However, 
because p$ falls off as a~ 3 immediately, compared to the a -4 fall-off for radiation, the result 
can be shown to be unacceptable over- dominance of the scalar field energy density today, 
so we will not consider this case further (this case corresponds to the lower portion of Fig. 
HD. The alternative, anharmonic limit, (j » 1 (second two terms dominant) yields 

r,] (1 + = a 4 Q, (3.3) 

with solutions 

f a^Q 1 / 3 « 1, for a « 1.5 x 10 A X 1 /% 1 , 
m « (3-4) 

1 2 1 / 4 «Q 1 / 4 » 1, for a » 1.5 x lO^V . 

The small a of the first case leads to a flatter potential (see Fig. P, with a minimum at a 

relatively large value of 0, but a small rj. The second case leads to a deeper potential, with a 

smaller, universal value at the minimum of (f>/M » (24/A) 1/4 (pv/6M 4 ) 1 / 4 , but a large value 

of r\. For these solutions, generally rji/(l + Pi) = rji/(l + 3a 2 ) « 1 for A and ty of order 

unity (see Eq. ( |3.6| ) below). 



B. Post-inflation evolution of scalar field 

At the end of inflation when reheating occurs, $ ~ 0, p — 3p = 0, and p$ ~ V(®i) = 
(^m 2 cj) 2 I + |A0 4 )/(1 + acfii/M) 2 . For simplicity, we assume that the vacuum energy density 
is totally converted to radiation energy, so that Pr = pv- Since we are considering the limit 
(i » 1 we can approximate Vi(§) ~ f^ 2 (^\(j) i ). Then using Eq. ( |3.3|) , we find that, right 
after inflation ends, 

Xi = (p^/r 2 p R )i « % 4 /2o 4 g « m /(2 + vi), (3.5a) 
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H 2 « ( Pv /QM 2 )(1 + + l - m y l . (3.5b) 



However, when the false vacuum energy density is converted to radiation and p — 3p — > 0, 
the minimum of the potential V goes to $ = 0. Whether $ begins to oscillate or rolls slowly 
down the potential V depends on whether the period of oscillation is short or long compared 
to the expansion time scale, where (period) 2 w ^(dV/d^)^ 1 ps 4>(d^/d4>) 2 (dV(^)/d4>)~ 1 . In 
this case, using d$/d(j) (Eq. (|2~?1)), we obtain (period) 2 ps m~ 2 (l + r]j)(l + /? + rji)(l + 2£j + 
fjiCi)' 1 - Since we have argued that r]i/(l+ ft) << 1, where ft = 3a 2 , and » 1, we then 
find the ratio 

_ (period) 2 A+ja 2 ^ / ^ \ 

^-F^r^A^i" (3 - 6) 

For the solution with rji « 1, 1Z « 1, so oscillations begin immediately, but with 
sufficiently small amplitude that we can approximate / m 1 right away. Then, initially, 
p$ ~ E oc a~ 4 , since the oscillations are anharmonic, and m oc a -1 . Since Cm oc m oc a -2 , 

1 /2 

the transition to harmonic oscillations (H) at C m ps 1 occurs at a# ~ a/. In the harmonic 
regime, p$ cx a -3 and m oc a~ 3//2 . It then follows that the present ratio of scalar to radiation 
energy densities is given by 

Xo = (p*//~ 2 Pfl)o ~ X/C/~ 1/2 («o/a/) ~ {2\y 1/2 (a/y)(a /a I ), (3.7) 



and the present amplitude of G-oscillations is r] ps* ^/C/^ 2 ( a // a o) 3 ^ 2 - Substituting Eq. (|2.21|) 
for ao/ai into Eq. ( |3.7| ), we obtain for the present scalar density parameter 



Q%h 2 = Q° R h 2 X o « 10 22 (2/A) 1 / 2 (a/ 2 /)(T v /10 14 GeV), (3.8) 

where we have used Eq. ( |2.20| ) for the radiation density parameter and have assumed 
T = 2.73° K. If G-oscillations are to account for the energy density needed to close the 
universe without other forms of dark matter, then we require Q% ~ 1. For dimensionless 
coupling A of order unity and Ty ps 10 14 GeV, the values of the scalar mass m = M/y and the 
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coupling constant a that will meet this constraint are shown in Fig. |^. If the universe is flat, 
then h < .65 is needed to satisfy age constraints from globular cluster ages. In Fig. 2, the 
solid dark line corresponds to Q%h 2 = 0.25, e.g., Q% = 1 and h = 0.5. The curve shows that 
Q% near closure density is possible for a wide range of parameters, including the range where 
dimensionless parameters, such as A and a of are of order unity. Notice in Fig. |2| that the 
present amplitudes of G-oscillation are extremely small, despite contributing a potentially 
critical energy density. Also plotted is the curve for Q%h 2 = 1 (dashed curve, parallel to 
dark curve), which is, roughly, the upper bound on Vt%h 2 consistent with observations. 

For the case rjj » 1, corresponding to the right-hand strip of Fig. the scalar energy 
density is comparable to the radiation energy density just following inflation (Eq. (|3.5a|) ), 
and the ratio 1Z (Eq. ( [3.5|) ) is of order unity. As a result, oscillations do not start imme- 
diately, instead, the scalar field rolls slowly down the potential, while the radiation energy 
density falls off rapidly. Eventually, 77 (and thus TV) becomes small enough that oscilla- 
tions commence, first anharmonically, and then harmonically as before. But because of 
the intervening decrease in the radiation energy density, the resulting scalar energy density 
overdo minates the universe today. This effect causes the curves of constant VL%h 2 in Fig. |2| 
to turn up as they approach the line where r]i ~ 1 (a ~ 10 4 A 1 / 4 ty 1 ). 



IV. QUADRATIC ALLY COUPLED OSCILLATING G 

A. Oscillations generated by inflation 

In the quadratically coupled case, the effective inflationary scalar potential (Eq. (|2.11| )) 
is given by 

n ' (l + ^ 2 /M 2 ) 2 1 ; 

Note that as — > 00, the potential tends to a constant, AM 4 /4£ 2 . Fig. [I] plots Vj($) 
schematically. The minimum of the potential occurs at 
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m = 2i 2 Qr { (4.2) 



where Q is given by Eq. ( gl9|) , and 



rB l-W>v/6Af]- >a (4 . 3) 

i - £y A 1 

If r < 0, the minimum occurs at rjj = 0, and no oscillations are generated by inflation. If 
r > 0, but both numerator and denominator in Eq. (|4.3| ) are negative, then r]i = 2£ 2 Qr 
is a maximum, while r)i = is a minimum, again resulting in no oscillations. These cases 
correspond to the region labelled "no oscillations" in Fig. |3|. For the quadratically coupled 
model, the anharmonicity parameter is given by ( = \\y 2 rj/^. Note that for oscillations 
to occur, the numerator of Eq. ( |4.3| ) must be positive, thus 24^y 2 (Hj/M) 2 = Q/r > 1, 
corresponding to the anharmonic regime. The larger (j is, the closer r is to unity. 

B. Post-inflation evolution of scalar field 

In this case, we find at the end of inflation, 

Xi « V 2 i/^ 2 Q « Vi, (4.4a) 

H 2 « (p v /QM 2 )(l + Vl r)(l + Vl )- 2 . (4.4b) 

The effective oscillation period is (period) 2 ~ m~ 2 (l + r]j)(l + j3i + r] I )(l + 2( I — r/j) -1 , Using 
Eq. (|4.2|) and the fact that Q » 1, we find that the ratio 1Z is 

^ l + fr + vi „ i + (igg + i) m 

24£r 24£ 1 ' ' 

There are three possible post-inflation evolutions. 

(a) Immediate Oscillations. Oscillations begin immediately when 1Z « 1, which cor- 
responds to 1/24 < £ < 2.9 x 10 7 A 1 / 2 . In this range, r/j « 1, so we approximate / ~ 1, 
and p$ << pr. However, although (/ >> 1, /?/ = 12£r)i ranges from 10~ 18 to 10 11 , with 
Pi « 1 for £ ^ 3 x 10 5 A 1 / 3 *v 

r 1 / 3 , so we must consider two separate cases (for the range 
15 



of parameters considered, Pi/Ci ~ 24£ 2 /y 2 << 1 ). Recall that f3 = |(1 + rj)/u, so that 
large j3 corresponds roughly to small u and to strong Brans-Dicke-like effects (we call this 
"strong non minimal coupling" or "strong coupling") and small (3 corresponds to large uj 
and to weak Brans-Dicke effects (we call this "weak coupling"). 

Case 1: f3i » 1: Because of the strong non- minimal coupling, the energy density falls 
off like matter rather than like radiation despite the anharmonic potential (Table |), so that 
initially E fa p$ oc a" 3 , while /3 m , £ m , Vm oc a~ 3 ^ 2 ; The transition to weak coupling (WC) 
((3 m fa 1) occurs at awe ~ ffl 3 ai, whereupon standard anharmonic oscillations lead to E oc 
a -4 and f3 m , ( m , rj m oc a~ 2 . Harmonic oscillations commence at a H fa Cwc a wc ~ Ci 2 ^]^ a i- 
During harmonic oscillations, E oc a -3 and f3 m , ( m , r) m oc a~ 3 . The result at the present 
epoch is 

Xo ~ XiCJ l,2 f3) ,2 {ao/ai) ~ Vi(7 1/2 ' ^{ao/aj), (4.6a) 
Vo ~ 77iC! /2 & 1/2 (a//a ) 3 , (4.6b) 

and 

fi|/i 2 w 4 x 10 6 A- 3/2 (e 3 /l/)(r y /10 14 GeV) 5 . (4.7) 

Note that, for larger values of £ within the range quoted above, namely 7 x 10 5 A 5 / 13 < 
£ < 2.9 x lC/A 1 / 2 the scalar energy density can actually surpass the radiation energy density 
during the strong-coupling regime. For the smaller values of £, the scalar density approaches 
but does not exceed the radiation during the strong-coupling regime; it is only when har- 
monic oscillations begin that the scalar energy surpasses the radiation. These two sub-cases 
are shown as curves (a) and (b) in Fig. |j. 

Case 2: Pj << 1: In this case, weak coupling holds right after inflation. Standard 
anharmonic oscillations begin immediately, E oc a -4 and j3 m , ( m , r] m oc a~ 2 , with harmonic 

1/2 

oscillations taking over when an ~ (j aj (curve (c) in Fig. [|). The evolution then proceeds 
as in the previous case, with the result 
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Xo ~ mCi 1/2 (a /a I ), (4.8a) 

Vo « ^(^'(az/ao) 3 , (4.8b) 

and 

ft°/i 2 « 1.6 x 10 14 A" 1 (e 3/2 /l/)(r y /10 14 GeV) 3 . (4.9) 

f&j Deferred oscillations. For £ << 1/24, r/j << 10~ 17 , xi « 1, but TZ » 1 and 
>> 1. This corresponds to small amplitudes for 0, though still in the anharmonic regime, 
but since m/H « 1 at this epoch, the scalar field first rolls slowly down the potential. Thus 
we can ignore <E> in Eq. ( |2.10b| ), and solve the approximate "slow-roll" evolution equation 

= (d(/)/d$)$ w -CiH)-\d(t)/d^fdV^)/d(t) « -(m 2 2 /3#)(l + 2(). (4.10) 

Radiation strongly dominates the energy density following inflation, so H ~ l/2i. Noting 
that C oc 2 , we rewrite Eq. ([4.1 0|) in the form £ « — |m 2 /;C(l + 2£). For mi ~ m/iJ << 1 



we find the approximate solution ( ta 0/(1 + |C/m 2 i 2 ). Noting that m 2 £/ ps 12£(py/6M 2 ) 
and defining r = t(py/6M 2 ) 1 / 2 , we obtain 

ft » (2t)- 1 (2m 2 C)" 1 w [1 + 16£r 2 ]/96£r 2 . (4.11) 

Thus oscillations begin when r w (SO^)" 1 / 2 . At this point, p# ~ 2p/r~ 2 ps 20£p/, p$ ~ 
A0 4 ps (25/36)p$, and a ps a//(20£) 1//4 . Following the oscillations from this point onward, 
we find 

n%h 2 Pd 10 13 A- 1 (^ 3/4 / 2 /)(T y /10 14 GeV) 3 . (4.12) 

This corresponds to curve (d) in Fig. f|. 

(c) Lingering inflation. For £ >> 3 x 10 8 A 1 / 2 , rji » 1, we have p$ ps (A/4£ 2 )M 4 , 
Xr ~ r/j >> 1, and 1Z ~ 7p/2 >> 1. In this case, the scalar-field energy density dominates, 
and the scalar field rolls slowly down the flat wing of the post-inflation potential (Figure p. 
In this case, the slow-roll evolution equation gives 
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<j> w -(3H)-\d<f>/dQ) 2 dV($)/d<i> « -i(A/6) 1/2 r 2 M 3 0-\ (4.13) 

o 

where we have used H pa (A/24) 1/,2 £ _1 M, and /3, 77, £ >> 1. The solution is 77/ — 77 « 
|(A/6) 1 ' 2 ^~ 1 M(t — £/). Oscillations in $ commence when 77 ~ 1 << r/i, i.e. when H(t — tj) pa 
I777. During this period, ~ constant, a oc e m , and the radiation density deflates by 
e -4H(t-ti) w e _6r?/ . This is a kind of "lingering inflation", driven by the approximately 
constant scalar energy density on the flat portion of the potential at large <fi (Figure p. 
Eventually oscillations begin and the evolution proceeds as in the previous case, except 
that the lingering inflation leads to an effective initial xi larger than before by the factor 
e 6??j >> fj^g resu it j s a much larger value for the scalar- field energy density: 

n%h 2 w 4 x 10 6 A- 3 / 2 (e 3 /7/)e( 4><10 ~ 16A ~ 1 « 2 ^)(T y /10 14 GeV) 5 . (4.14) 

Because ^ > 3 x 10 8 A 1 / 2 in this case, we see that the scalar field overdominates the present 
energy density by many orders of magnitude. This range of values of £ is therefore unviable. 

As discussed for the linear case, Vt%h 2 ~ 0.25 corresponds to the case where the oscillation 
energy in $ is near the closure density, assuming no additional dark matter. Imposing 
Vt%h 2 pa 0.25 and assuming A and Ty/10 14 GeV of order unity, we obtain from Eqs. ( (4.7|) , 
(4.9), ( f4.12|) , and ( 4.14Q the constraint on the scalar mass m = M/y and the coupling 



constant £ shown in Fig. ^. As with the linear case, we find that Q$ near closure density is 
possible for a wide range of parameters, including the range where dimensionless parameters, 
such as A and £ of are of order unity. Also shown is the corresponding curve for VL%h 2 = 
1, roughly the largest value compatible with observations.. Not surprisingly, given the 
quadratic nature of the non-minimal coupling, the amplitude of G-oscillations at the current 
epoch is extraordinarily small, despite their ability to generate a closure energy density. 

V. CONSTRAINTS FROM LOCAL EXPERIMENTS 



We turn now to a discussion of bounds on these models imposed by local laboratory 
and solar-system experiments. Because of the presence of the potential V($), there is a 
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length scale i over which modifications to general relativity may be important. For scales 
much greater than £, the theory is equivalent to general relativity. Numerous tests of general 
relativity and of the inverse-square law of gravity using a variety of techniques have been 
performed which cover a range of lengths from solar-system scales down to centimeters. In 
this section we present a general method for using the results of such experiments to bound 
the parameters of oscillating G theories. We will find that the most viable models from a 
cosmological viewpoint are comfortably consistent with all local experimental limits. 

Comparison with experiment is most easily done in the physical representation. We 
expand the physical metric and the scalar field about their asymptotic values, = rj^+h^, 
and = 4>o(t) + f, and write /(</>) = f + f' ip + ^f^ip 2 , and substitute into Eqs. Q2.4|), 



keeping terms to first order in and ip. In eliminating the zeroth order terms in the 
field equation for 0, corresponding to the cosmological solution for O > we note that the 
coordinate t used here is related to the time t^w of the physical Robertson- Walker metric 
by g^rw = (1 — ^h 00 )dt. The result is 

- , + (V - = 4 ydW'" " 3P) + S*' + 5*"** (5 ' 1) 



where 



~ 2 _ m 2 f 2 />o(/o-/W(l + 6M 2 /;) 



4S^W^-^ + 3MVM}. (5.2) 

We are interested in approximately static solutions corresponding to a gravitating mass 
such as the Sun or Earth, approximately at rest (v m 0,p << p). Thus we can drop the 
terms involving time derivatives. Using Eq. ( |2.3| ), and averaging the oscillations of 0o over 
several periods, we find 

whose solution is 

ip = 4M 2 f^e-™, (5.4) 
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where p = / ~ 1 GM source is the effective gravitational mass of the gravitating body. 

Using the approximations /. o ~ 0, □/ ~ /gV 2 v9, T 00 ~ p, Tij ~ 0, and working in a 
gauge in which hf ^ — . = (/o//o)</>,i5 we obtain the field equations 

V 2 (/i 00 - UUfM = ~p/2M 2 f - ( P(f> + 3^)/2M 2 / , (5.5a) 

V 2 (^- + UUh)Sij(p) = S ijP /2M 2 f - Syfo - ^)/2M 2 / 0l (5.5b) 

where p^ = + ^(0o) and p<f, = |0q — V(<j>a). Dropping the p^ and terms for the 
moment, we obtain the solutions 

ft »-T( 1 + 3T^)- (5 - &) 

h -> = t s " i 1 ~ srs?"*) ■ (5 - 6b) 

The effective mass m leads to a range vrC x for modifications of Einstein gravity. The 
"Newtonian" potential /loo can be written in the more familiar form /loo = 2G e ///x/r, where 
G e fj has the effective form 

G * = ( 1 + 3T^^')- (57) 

Because of the Yukawa-type correction to the potential, for r << rh~ x , G e ff ~ (4 + 2a>)/(3 + 
2cu), while for r >> m -1 , G e ff ~ 1. Thus, for r >> mT 1 , the theory is equivalent to general 
relativity. Experimental tests of classical general relativity can in principle be used to place 
bounds on m and on the coupling parameters a and £ of the two oscillating G models. 

A. Light deflection and time delay 

To lowest order in /i/r, the physical metric obtained from Eqs. (|5.6| ) can be written in 
the form 

2u w 2u e-™\ 

20 



+ (5.8b) 

Notice that, in the line element, ds 2 = g IM/ dx lx dx v , the exponential terms can be extracted 
as an overall conformal factor. Hence, there is no direct effect on null geodesies, and the 
deflection and time delay of a signal is proportional to /x, with no explicit dependence on m. 
However, \i represents the mass as measured at spatial infinity using Kepler's law. In solar 
system applications, such as the deflection of light, the solar mass actually used is measured 
by using Kepler's law applied to the Earth's orbit. If the range fh~ l is small compared to 
the Earth's orbit, i.e. m~ l « 10 8 km, corresponding to in » 10 -27 GeV, or an oscillation 
frequency v » 1CT 2 Hz, the inferred solar mass M & is the same as /i in Eq. ( |5.8|) , and 
consequently, the predicted time delay and light deflection in terms of the measured solar 
mass are indistinguishable from Einstein relativity. This corresponds to y « 10 45 , which 
easily includes the masses of interest for cosmology. Conversely, if fhr 1 » 10 8 km, the 
inferred mass is M & = /i(4 + 2u)/(3 + 2u), and the prediction for time-delay and light 
deflection approaches the Brans-Dicke form, proportional to (3 + 2u)/(A + 2a;) . 



B. Orbits 



Using Earth-LAGEOS-Lunar data which compares GM ffi inferred from Earth-surface 
gravity measurements to the values inferred from the study of the orbits of the LAGEOS 
satellites and the Moon, tight bounds can be obtained on the variation in G in the 10 4 - 



10 m range. From a rough fit to the appropriate portion of Fig. 1 in [11], we find that the 



coefficient in the exponential correction in Eq. ( |5.7|) is constrained to be 

(3 + 2uj)- 1 < 0.03(z//MHz) L1 , (5.9) 

for frequencies in the 10 kHz to 1 Hz range, corresponding to values of y between 10 38 and 
10 42 . For the higher frequencies of interest to our cosmo logical models, the bound is easily 
satisfied. 
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C. Tower and laboratory measurements 



Tower measurements of the variation in G place the limit (3 + 2u)~ 1 < 2 x at 
a range of 100 meters i.e. v = 3 MHz, or y = 10 36 . Stronger limits are obtained from 
laboratory tests. From a null test of Gauss's law at 1.5 meters using superconducting 
gravity gradiometers [O (y = 200 MHz, y = 10 34 ), we find the bound 

(3 + 2U)- 1 < 3 x 10~ 5 £ 2 e 3/i , (5.10) 

where I is the range in meters. From a test of the inverse square law at 2 centimeters using 
a torsion balance |l3j {y = 15 GHz, y = 10 32 ), we obtain 



(3 + 2u)- x < 3.4 x 10~ 2 £ 2 e 0M/e . (5.11) 

The bounds relax exponentially with increasing mass and frequency (decreasing range), 
and quadratically with decreasing mass and frequency (increasing range). Notice that the 
conversion from range I to the variable y is £/lmeter = 1.14 x 10~ u y(m/m). Since the range 
of values of y of interest to cosmology is y < 10 30 , only the bound from the shortest-range 
experiment, Eq. ( |5.11|) will be meaningful. 



D. Application to specific models 

Linear Coupling. With small amplitude oscillations at the present epoch (r] « 1), we 
find from Eqs. Q and ( |OD 

uj = l/2a 2 , m 2 = m 2 /(l + 3a 2 ). (5.12) 

For y < 10 30 , Eq. (|5.11|) yields the bound ay < 4 x 10 31 . This bound appears in the upper 
right corner of Fig. 

Quadratic Coupling. In this case, to first order in r/ and /3 , we have 

uj = 3/2p , m 2 = m 2 (l - 5r/ - 2(3 ). (5.13) 
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But because r] and (3q are so small (see Fig. |3|), u> is sufficiently large to be compatible with 
any experimental bounds; in this case, despite the massive scalar field, general relativity is 
a strong attractor at the present epoch. 



E. Effective cosmological constant 

We now consider the terms in Eqs. (5.5) that involve the energy density and pressure 



of the oscillating scalar field. For harmonic oscillations, averaged over a period, (p^) ~ 
constant, varying only on a Hubble timescale, while (p^) ~ 0. The resulting spatially uniform 
cosmological scalar energy density acts like an effective cosmological constant, resulting in 
contributions to the local physical metric of the form 

gm ~ _i _ Qj# o V , (5.14a) 



0y « - fijtffr 2 ) . (5.14b) 

These terms will modify the gravitational dynamics of local systems, ranging from the solar 
system to clusters of galaxies. In the case of the usual cosmological constant A, observa- 
tional bounds from such systems are generally not sufficient to constraint the energy density 
represented by A to values below the critical energy density ||. Thus those observations are 
consistent with our desire to have Q% ~ 1. 



VI. CONCLUSIONS 

The scenario we have investigated in this paper links together several plausible compo- 
nents. First, we have assumed some modification of standard, Einstein gravity in the form 
of a non-minimally coupled scalar field. Modifications are required in almost any attempt to 
obtain a unified quantum theory of particle physics and general relativity, including super- 
gravity and superstring theories. The specific non-minimal coupling considered in this paper 
is similar to the coupling of the dilaton. Second, we have assumed that the non-minimally 
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coupled field is massive and has typical self-interactions of quantum scalar fields. Third, we 
have assumed inflation, which is well-motivated for many independent reasons. 

We have shown that, taken together, these components lead naturally to a scenario 
in which G oscillates about its mean value. Inflation shifts the scalar field away from its 
minimum. At some time after inflation, the field begins to oscillate. If the mass of the scalar 
field is typical of elementary particles, m > 1 GeV or so, the oscillations are high frequency 
compared to the expansion rate of the Universe. The amplitude is exponentially small, 
though, since the Hubble expansion steadily redshifts the kinetic energy density between 
inflation and the present epoch. For a wide range of model parameters, the combination 
of the high frequency and small amplitude that typically arises passes all known tests from 
cosmology and general relativity. 

Our analysis applies to models in which the $ field evolves purely according to the 
classical equations of motion ( |2.10b| ) and the oscillation energy dissipates through Hubble 
red shift. Since the $-field is nonminimally coupled, it is possible in many theories for 
the $-field energy to dissipate also through decay into other particles via gravitational 
interactions. Decay commences once the age of the Universe exceeds the lifetime of the 
$-field. The classical calculations leading to Figs. 2 and 3 significantly overestimate the 
oscillation amplitude and energy today. The possibilities depend upon the details of the 
model and lie beyond the scope of this paper. 

The provocative possibility raised by the analysis is that the oscillation energy in G could 
account for a large fraction of the energy density needed to close the Universe. Here, even 
the most extreme case, where oscillation energy is the only nonluminous energy density, 
has been shown to be viable for a wide range of parameters, including the range for which 
dimensionless parameters are of order unity. A future challenge is to determine if oscillation 
energy by itself or combined with conventional dark matter can account for large-scale struc- 
ture formation. Perhaps these investigations will motivated new experimental approaches 
for measuring high frequency oscillations of Newton's constant, or other local effects of such 
theories. 
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FIGURES 



FIG. 1. Schematic plot of effective inflationary scalar potential during and after inflation 
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FIG. 2. Allowed values of y = M/m and a in the linearly coupled model. Hatched line 
corresponds to bounds from laboratory inverse square-law tests. The solid dark line corresponds to 
n%h 2 = 0.25, the case of a critical density of scalar field energy today and h = 0.5. Also shown is 
VL%h? > 1 (dashed line and shading), roughly the upper limit consistent with present observations. 
Dotted lines are lines of constant amplitude of present-day G-oscillations (770). Also shown (dashed 
lines with shading) are the excluded regions in which ?/; > 1 or (/ < 1. 
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FIG. 3. Allowed values of y = M/m and £ in the quadratically coupled model. Here, the 
oscillation amplitude (770) is so tiny that laboratory experiments provide no useful bounds. The 
excluded regions are where the field is not displaced from the minimum of the effective potential 
during inflation (bottom shaded region) or where there is lingering inflation due to the $ field such 
that oscillation energy ultimately overdominates the universe. 
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FIG. 4. Schematic evolution of scalar energy density, pf 2 , in quadratically coupled model for 
various ranges of coupling parameter £ and mass m, all leading to the same scalar energy density 
today. Superimposed are solid curves showing the evolution of the baryon density, ps, and the 
radiation energy density, pR. The scalar energy density corresponds to the various dashed curves: 
Curves (a) and (b): Oscillations begin immediately after inflation, but in the strong nonminimally 
coupled regime, where pf 2 oc a~ 3 . Depending on value of £ scalar energy density may grow to 
exceed that of the radiation density, e.g., curve (a), before entering the weakly-coupled regime 
where pf 2 begins to evolve as a -4 . Curve (c): Immediate oscillations in weakly-coupled regime 
fall off as radiation from the start. Curve (d): Deferred oscillations cause scalar energy density to 
remain roughly constant initially, then to fall off like radiation. Eventually harmonic oscillations 
take over in all cases, leading to matter-like falloff. In all cases, the late, weakly-coupled behavior 
switches to harmonic oscillations with pf 2 oc a~ 3 . 
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TABLES 



TABLE I. Values of 7 and 5 for scalar-field oscillations varying according to E oc a 37 and 

(p m OC a~ 3S 





Linear 


Quadratic 








/3m « 1 


An » 1 


Values of 7 


Harmonic Cm « 1 


1 


1 


2/3 


Anharmonic Cm » 1 


4/3 


4/3 


1 


Values of 5 


Harmonic Cm << 1 


1/2 


1/2 


1/3 


Anharmonic Cm >> 1 


1/3 


1/3 


1/4 
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